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ABSTRACT 

Dark energy dominates the energy density of our Universe, yet we know very little about its nature and origin. 
Although strong evidence in support of dark energy is provided by the cosmic microwave background, the relic 
radiation of the Big Bang, in conjunction with either observations of supernovae or of the large scale structure of 
the Universe, the verification of dark energy by independent physical phenomena is of considerable interest. We 
review works that, through a wavelet analysis on the sphere, independently verify the existence of dark energy by 
detecting the integrated Sachs- Wolfe effect. The effectiveness of a wavelet analysis on the sphere is demonstrated 
by the highly statistically significant detections of dark energy that are made. Moreover, the detection is used 
to constrain properties of dark energy. A coherent picture of dark energy is obtained, adding further support to 
the now well established cosmological concordance model that describes our Universe. 
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1. INTRODUCTION 

Dark energy represents the energy density of empty space, providing space with an intrinsic nature to expand. 
The revolutionary idea that empty space contains energy was first proposed by Einstein, when he introduced 
the cosmological constant A to reconcile his general relativistic model of the Universe with his prejudice towards 
a static universe. Observational evidence^ later contradicted Einstein's assumption, instead showing that the 
Universe is expanding. Einstein subsequently referred to his introduction of the cosmological constant, which 
precluded him from predicting the expansion of the Universe, as his "greatest blunder" . Of course, only Einstein's 
greatest blunder could prove to be a remarkable insight which is today supported by strong observational evidence 
(although it was postulated for the wrong reasons initially). We review in this paper a method that has been 
used to confirm independently the existence of dark energy and to constrain its properties. 

The integrated Sachs- Wolfe (ISW) effect^ provides an independent physical phenomenon that may be used 
to detect and probe dark energy. This phenomenon predicts a secondary contribution to the anisotropics of the 
cosmic microwave background (CMB), the relic radiation of the Big Bang. CMB photons undergo an energy 
shift due to the evolution of gravitational potentials as they travel towards us from emission shortly after the 
Big Bang. This physical phenomenon is known as the ISW effect. The secondary contribution to the CMB due 
to the ISW effect is only present if the Universe deviates from matter domination due to either spatial curvature 
or dark energy. Strong constraints have now been placed on the fiatness of the Universe,^ hence a detection of 
the ISW effect may be inferred as direct confirmation of the existence of dark energy. It is difficult to isolate 
the ISW contribution to CMB anisotropics directly. Instead, as first proposed by Ref. 4, it is possible to detect 
the effect by searching for a correlation between the CMB anisotropics and the large scale structure (LSS) of 
the Universe. A positive large-scale correlation will be induced by the ISW effect as a consequence of decaying 
gravitational potentials due to dark energy. First attempts to detect the ISW effect using Cosmic Background 
Explorer-Differential Microwave Radiometer (COBE-DMR) CMB observations failed, concluding that greater 
sensitivity and resolution than that provided by COBE were required.^ Fortunately the Wilkinson Microwave 
Anisotropy Probe (WMAP) mission soon provided suitable CMB data.^'^ Correlations indicative of the ISW 
effect were first detected by Ref. 8 using the WMAP data and the NRAO VLA Sky Survey (NVSS)^ radio galaxy 
survey, and have since been detected between various releases of the WMAP data and a large number of tracers 
of the LSS of the Universe. 
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Many different analysis techniques have been employed to detect the ISW effect, each of which has its own 
merits and limitations. We focus here on wavelet-based techniques. Wavelets provide an ideal analysis tool to 
search for the ISW effect. This is due to the localised nature of the physical phenomenon, in both scale and 
position on the sky, and the simultaneous scale and position localisation afforded by a wavelet analysis. Since 
detections of the ISW effect are limited by the proportion of the sky observed, it is desirable to have as great a 
sky coverage as possible. In this near full-sky setting the geometry of the celestial sphere on which observations 
are made should be taken into account. Consequently, wavelet analyses on the sphere are required. A variety 
of different methodologies have been proposed to construct a wavelet formalism on the sphere (see e.g. Ref. 10 
for a review). Attempts to detect the ISW effect have predominately applied the methodlogy proposed by 
Ref. 100 The use of wavelets on the sphere to search for correlations indicative of the ISW effect was pioneered 
by Ref. 12, using the axisymmetric spherical Mexican hat wavelet. Correlated features induced by the ISW 
effect may not necessarily be rotationally invariant; indeed, it is known that statistically isotropic Gaussian 
random fields are characterised by local features that are not rotationally invariant. The analysis performed 
by Ref. 12 was extended to directional wavelets by Ref. 14 in order to probe non-rotationally invariant features. 
A fundamentally different analysis using wavelets on the sphere has been performed by Ref. 15 recently. This 
method employs steerable wavelets on the sphere to probe the morphology of local features. All of these works 
examined the WMAP and NVSS data and made significant detections of the ISW effect, thereby confirming the 
existence of dark energy. The effectiveness of the wavelet approach is highlighted in these analyses by the highly 
statistically significant detections of the ISW effect that have been made. Moreover, in the first two of these 
works (Refs. 12, 14) the detection of the ISW effect was used to constrain dark energy properties. In this paper 
we review the aforementioned works that employ a wavelet analysis on the sphere to detect and constrain dark 
energy. 

The remainder of this paper is organised as follows. In Sec. [2] we present the cosmological foundation on 
which the reviewed works are based. The cosmological model that describes our Universe is discussed briefly, 
before the origin of the ISW effect is described. In Sec. [3] the wavelet methodology on the sphere applied to 
detect the ISW effect is presented. Methods to detect the ISW are outlined in Sec. [H and positive detections of 
the effect are presented. Dark energy properties are also constrained in Sec. [H Concluding remarks are made in 
Sec. El 

2. COSMOLOGY 

Modern cosmology has converged on a concordance model only recently. The so-called A Cold Dark Matter 
(ACDM) cosmological concordance model describes many independent observations of our Universe very well. 
For example, the model describes the relic radiation of the Big Bang {i.e. the CMB), the LSS of our Universe and 
the accelerating expansion of our Universe apparent from observations of supernovae. We describe this model 
very briefly here, before using it to explain the origin and observational implications of the ISW effect. 

2.1 Cosmological concordance model 

The ACDM model is characterised by a Universe consisting of ordinary baryonic matter, cold dark matter and 
dark energy. Current estimates place the relative contributions of these components at 4%, 22% and 74% of the 
energy density of the Universe respectively.^ Cold dark matter consists of non-relativistic, non-baryonic patricles 
that interact gravitationally only. Although dark matter has yet to be observed directly, its presence has been 
inferred from a large range of observations, such as galaxy rotation curves. Dark energy represents the energy 
density of empty space and may be modelled by a cosmological fluid with negative pressure, acting as a repulsive 
force counteracting the attractive gravitational nature of matter. In general, the dark energy fluid with pressure 
p and density p has an equation of state p = wp, where w = —1 corresponds to the case of a cosmological 
constant. Evidence for dark energy is provided directly from observations of the statistical properties of the 
CMB anisotropics, together with either observations of Type la Supernovae^^' or of the LSS of the Universe 
{e.g. Ref. 19). Although strong evidence in support of dark energy exists, a consistent model in the framework 
of particle physics is lacking. Indeed, attempts to predict a cosmological constant using quantum field theory 

*A new wavelet construction on the sphere called needlets has also been used recently to detect the ISW effect,^^ 
although we do not focus on this work in any further detail here. 



of the vacuum energy density arising from zero-point fluctuations predict a value that is too large by a factor 
of ~ 10^^^. Dark energy dominates the energy density of our Universe and yet we know very little about its 
origin and nature. The verification of dark energy by independent physical methods is therefore of considerable 
interest, hence the focus here on detecting dark energy through the ISW effect. Moreover, independent methods 
may provide more sensitive probes of certain properties of dark energy. 

Now that the fundamental components of our Universe have been described, we turn our attention to a 
dynamical view of the Universe. Astronomical observations made by Hubble^ (and since more accurately by 
others) showed that our Universe is expanding. Moreover, distant regions of our Universe are receding from 
us at a rate proportional to their distance. By reversing time and extrapolating into the past one recovers the 
Big Bang model, consisting of a universe that began from a singularity and expanded and evolved into the 
Universe observed today. However, the Big Bang model as stated currently has a number of problems. One 
of these problems is that the model contains no mechanism to generate the density perturbations required to 
explain the structure of the Universe that is observed today, such as galaxies and clusters of galaxies. The 
theory of infiation provides a mechanism to solve this (and other) problem(s) of the original Big Bang model. 
Inflation proposes that at a very early stage the Universe underwent a phase of rapid exponential expansion, 
increasing in size by many orders of magnitude. The precise details of inflation are an active area of research, 
with many competing inflationary scenarios. In essence, a period of rapid expansion could occur if the Universe 
underwent a phase-transition, perhaps associated with symmetry breaking of the grand unifled fleld, generating 
a false vacuum with a vacuum energy density acting like a cosmological repulsion. The details of a scalar fleld 
with negative pressure which may drive such an expansion have yet to be identifled. The inflationary scenario 
provides a mechanism by which initial perturbations in the primordial fluid may be generated, thereby providing 
the seeds of structure formation. Quantum fluctuations in the very early Universe, or more precisely in the scalar 
fleld driving inflation, would be blown up to macroscopic scales by inflation, essentially freezing the quantum 
fluctuations into the Universe. The resulting not quite perfectly homogeneous state may then evolve under 
gravitational instability to produce the structure of the Universe observed today. 

With the components and dynamics of the Universe now discussed, we are in a position to consider the flnal 
piece of the jigsaw required to understand the ISW effect: the CMB. The temperature of the early Universe 
was sufficiently hot that photons had enough energy to ionise hydrogen. Compton scattering tightly coupled 
photons to electrons, which were in turn coupled to baryons through electromagnetic interactions. Scattering 
happened frequently and hence the mean free path of photons was extremely small. At this time the Universe 
therefore consisted of an opaque photon-baryon ffuid. As the Universe expanded, this primordial plasma cooled 
until the majority of photons no longer had sufficient energy to ionise hydrogen. Consequently, the fraction 
of ionised electrons dropped, the photons decoupled from the baryons, and the Universe became essentially 
transparent to radiation. This process, know as recombination^ occurred when the temperature of the Universe 
fell to approximately 3000K, corresponding to time of approximately 400,000 years after the Big Bang (a small 
period relative to the age of the Universe, estimated at 13.7 billion years). The photons were then free to 
propagate largely unhindered through space and may be observed today as the CMB radiation. The CMB is 
highly isotropic, with anisotropics originating from primordial perturbations at a level of 10~^ only. The CMB 
was first observed by Penzias and Wilson in 1965^^ and since by many other ground and satellite based telescopes. 
In Fig. [T] we show the temperature anisotropics of the CMB observered on the celestial sphere by WMAP. 

2.2 Integrated Sachs- Wolfe effect 

As CMB photons travel towards us from emission shortly after the Big Bang, they pass through gravitational 
potential wells due to the LSS of the Universe. The photons suffer blue and red shifted as they fall into and 
out of potential wells respectively. If the gravitational potential evolves during the photon propagation, then the 
blue and red shifts do not cancel exactly and the photon undergoes an energy shift that produces a secondary 
contribution to the temperature anisotropy of the observed CMB. The ISW effect is the integrated sum of this 
energy shift along the photon path. Any recent acceleration of the expansion of the Universe due to dark energy 
will cause local gravitational potentials to decay. CMB photons passing through over dense regions of decaying 
potential suffer blue shifts, resulting in a positive correlation between the induced anisotropy and the local matter 



Figure 1. Full-sky map of the CMB anisotropies measured by WMAP (courtesy of the WMAP science team: 
[http: //map. gsfc .nasa.gov/ 1 ). The Mollweide projection is used here and subsequently to display data defined on the 
sphere. 



distribution. It can be shown^ that the temperature fluctuation induced in the CMB is given by 

A^H = ^|M = _2 J dvHv,u^), (1) 

where AT is the induced temperature perturbation, Tq is the mean temperature of the CMB, rj is conformal 
time, ^ is the gravitational potential and the dot represents a derivative with respect to conformal time. A 
point LO on the sky is represented in spherical coordinates as a; = (6>, with co-latitude and longitude cp. The 
integral is computed over the photon path from emission to observation. In a matter-dominated universe the 
gravitational potential remains constant with respect to conformal time {i.e. i> = 0), thus there is no ISW effect 
(the expansion of the Universe exactly offsets the gravitational infall of over-dense regions). However, if the 
Universe deviates from matter domination due to curvature or dark energy then an ISW effect is induced. As 
mentioned previously, strong constraints have been placed on the flatness of the Universe, thus any ISW effect 
is a direct indication of dark energy. 

It is difficult to isolate the ISW contribution to the CMB anisotropies directly. Instead, the effect may be 
detected by searching for a correlation between the CMB and the LSS of the Universe.^ The galaxy count 
distribution projected on to the sky 6^{z^u) may be used as a tracer of the local matter distribution S{z^uj). It 
is assumed that two relative fluctuations are related by the linear bias factor 6, i.e. 5^{z,u) = b{z) 6{z^ij), where 
z is redshift and the bias is assumed to be redshift independent. The galaxy source count fluctuation observed 
on the sky is therefore given by 

ANH=6 I dz^S{z,u;), (2) 

where dN/ dz is the mean number of sources per steradian at redshift z and the integral is performed from today 
to the epoch of recombination. 

It is possible to derive the theoretical correlation between the CMB temperature and galaxy count pertur- 
bations specified by Eqn. ([T]) and Eqn. ([2]) respectively. The correlation is represented in harmonic space by the 
cross-power spectrum Cf^ , defined by the ensemble average of the product of the spherical harmonic coefficients 
of the two signals observed on the sky: 

{AZAj,:,)=5ee'Snn'Cf^ , (3) 

where A^^ = (Y^^IA) are the spherical harmonic coefficients of A(a;), (•!•) denotes the inner product on the 
sphere, Y^n are the spherical harmonic functions for multipole ^gN, nGZ, \n\ < i and Sij is the Kronecker 
delta symbol. In writing the cross-correlation in this manner we implicitly assume that the galaxy density and 
CMB random fields on the sphere are homogeneous and isotropic. Representing the gravitational potential and 
the matter density perturbation in Fourier space and substituting Eqn. ([T]) and Eqn. ([2]) into Eqn. (jS]), it is 
possible to show {e.g. Ref. 21) that 

Cr = I2n Ho' J ^ Ajik) Frik) Fj{k) , (4) 



where Qrn is the matter density, Hq is the Hubble parameter, A^(/c) = Ps{k)/ 27r'^ is the logarithmic matter 
power spectrum, Ps{k) = {\5{k)\'^) is the matter power spectrum and the filter functions for the galaxy density 
and CMB are given by 

Ff(fc)=6 J dz^D{z)Mkr,{z)] (5) 

and 

Fj{k) = j dz^J,[kv{z)] (6) 

respectively. The integration required to compute Fj{k) is performed over z from zero to the epoch of recom- 
bination, whereas, in practice, the integration range for F^{k) is defined by the source redshift distribution 
dN/ dz. D{z) is the linear growth factor for the matter distribution: S{z^k) = D{z)S{k)^ with S{k) = S{0^k). 
The function g{z) = z)D{z) is the linear growth suppression factor and is the spherical Bessel function. 
We have represented in harmonic space the expected correlation between the galaxy source count and tempera- 
ture fiuctuations induced by the ISW effect. We next turn our attention to wavelets on the sphere as a potential 
tool for detecting this predicted correlation. 

3. WAVELETS ON THE SPHERE 

To perform a wavelet analysis of full-sky maps defined on the celestial sphere, Euclidean wavelet analysis must 
be extended to spherical geometry. The ISW analyses that we review adopt the wavelet transform constructed 
by Ref. 10. This wavelet construction on the sphere is derived entirely from group theoretic principles, however 
the formalism has recently been reintroduced in an equivalent, practical and self-consistent approach that is 
independent of the original group theoretic framework. We adopt this latter approach, with the extension to 
anisotropic dilations derived in Ref. 2313 In order to perform the wavelet analysis on high resolution cosmological 
data sets it is essential to apply fast algorithms to compute the analysis. The reviewed works make use of the 
fast algorithm derived by Ref. 23, which is in turn based on the fast spherical convolution derived by Ref. 24. 
We present here a brief overview of the wavelet transform on the sphere but refer the reader to Ref. 25 for a 
more detailed review of this wavelet methodology and corresponding fast algorithms. 

3.1 Wavelet transform 

The wavelet formalism on the sphere (S^) is constructed in an analogous manner to the construction of Euclidean 
wavelets. Affine transformations, such as translations and dilations, must be defined on the sphere. A wavelet 
basis on the sphere may then be constructed from affine transformations of a mother wavelet defined on the 
sphere. The wavelet transform is given by the projection of a function on to this basis, where the mother wavelet 
must satisfy an admissibility relation to ensure that perfect synthesis of the original function from its wavelet 
coefficients is possible. In the most general setting, the projection is given by a directional spherical convolution. 
For a given scale, the wavelet coefficients therefore live in the rotation group S0(3). 

In order to construct the wavelet framework on the sphere following the methodology outline above, it is 
necessary to extend the affine transformations of translation and dilation to the sphere. The natural extension 
of translations to the sphere are rotations. These are characterised by the elements of the rotation group SO (3), 
which we parameterise in terms of the three Euler angles p = (q^,/3,7)I1 The rotation of a function / on the 
sphere is defined by 

mfjiiu) = fip-'io), peS0(3). (7) 

The extension of the dilation operator to the sphere is not so simple and various operators have been proposed. 
Refs. 10,22 define a dilation operator that may be conceptualised through an association with the tangent plane 
at the north pole of the sphere, using the stereographic projection to map a function defined on the sphere to 

^Although the ability to perform anisotropic dilations is of practical use, we do not achieve a wavelet basis in this 
setting since perfect synthesis of the original signal is not possible. In any case, reconstruction is not required for this 
application. 

*We adopt the zyz Euler convention corresponding to the rotation of a physical body in a fixed co-ordinate system 
about the y and z axes by 7, /3 and a respectively. 



the tangent plane. The stereographic projection is defined by projecting a point on the sphere to a point on the 
tangent plane at the north pole, by casting a ray through the point and the south pole. The point on the sphere 
is mapped on to the intersection of this ray and the tangent plane (see Fig. [2j). We denote the stereographic 
projection operator by 11, with inverse 11"^, where these operators include the appropriate normalisation term to 
ensure that the L^-norm of functions on the sphere and tangent plane is preserved. In this formalism dilations on 
the sphere are constructed by first lifting the sphere to the plane by the stereographic projection, followed by the 
usual Euclidean dilation in the plane, before re-projecting back on to the sphere. We adopt the generalisation 
to anisotropic dilations on the sphere. The spherical dilation operator P(a, b) is defined as the conjugation by 
n of the anisotropic Euclidean dilation d{a, b) on the tangent plane at the north pole: 

V{a,b) =U-^ d{a,b)U, (8) 

where the Euclidean dilation operator cZ(a, b) has the appropriate normalisation factor to ensure that the L^- 
norm is preserved. Consequently, the L^-norm of functions is preserved by the spherical dilation as both the 
stereographic projection and Euclidean dilation operators preserve the norm of functions. The dilation of a 
function on the sphere / G L^(S^) is given by 

[D{a,b)f]{u;) = [X{a,b,e,^)]'/^ /K/a.i/b) , (9) 

where uja,b = {Oa,b^^a,b)^ tan(6>a,6/2) = tan(^/2)\/a^ cos^ (p -\- b'^ sin^ (p and tdin{(pa,b) = f tan((/?). For the case 
where a = b the anisotropic dilation reduces to the usual isotropic case defined by tan(6>a/2) = atan(6>/2) and 
(fa = (p. The A(a, 6, ^, cp) cocycle term follows from the factors introduced in the stereographic projection of 
functions to preserve the L^-norm. The cocycle of an anisotropic spherical dilation is defined by 

A(a,6,^,(^)^-- --—2 , (10) 

(A_ cos6> + A+) 

where 

A± = o?b^ ± o? sin^ ± 6^ cos^ (p . 

For the case where a = b the anisotropic cocycle reduces to the usual isotropic cocycle. Although the ability 
to perform anisotropic dilations is of practical use, one does not achieve a wavelet basis in this setting. In 
the anisotropic setting a bounded admissibility integral cannot be determined (even in the plane), thus the 
synthesis of a signal from its coefficients cannot be performed. Nevertheless, this is not of concern for the 
ISW analyses discussed herein since synthesis is not required. The projection of a signal on to basis functions 
undergoing anisotropic dilations may be performed in an analogous manner to the following discussion of the 
wavelet transform. However, since these basis functions are not wavelets we restrict the following discussion to 
isotropic dilations. 

A wavelet basis on the sphere may now be constructed from rotations and isotropic dilations (where a = b) 
of a mother spherical wavelet G L^(S^). The corresponding wavelet family {^pa,p = Tl{p)V{a^a)il) : p G S0(3), 
a G Mj} provides an over-complete set of functions in L^(S^). The wavelet transform on the sphere of a square 
integrable function / G L^(S^) is given by the projection on to each wavelet basis function in the usual manner: 

where dQ{Lj) = sin^ d^ d(p is the usual invariant measure on the sphere. The transform is general in the sense that 
all orientations in the rotation group SO (3) are considered, thus directional structure is naturally incorporated. 
It is important to note, however, that only local directions make any sense on S^. There is no global way of 
defining directions on the spherqj - there will always be some singular point where the definition fails. Although 
it is not of relevance to the ISW analyses reviewed herein, for completeness we also state the synthesis of a signal 
on the sphere from its wavelet coefficients: 

= I Mp) r ^ <(«,p) [n{p)LM{^) , (12) 

^SO(3) Jo <^ 



^ There is no differentiable vector field of constant norm on the sphere and hence no global way of defining directions. 



South pole 

Figure 2. Stereographic projection of the sphere on to the plane. 



where dg{p) = sin /3 da d/3 dj is the invariant measure on the rotation group SO (3) and the operator is defined 
by the action (L^g)^^ = gin/C^ on the spherical harmonic coefficients of functions g G L^(S^). In order to ensure 
the perfect reconstruction of a signal synthesised from its wavelet coefficients, the admissibility condition 



must hold for all £ G N, where {ipa)in spherical harmonic coefficients of V^a-^^ 

3.2 Correspondence principle 

The correspondence principle between spherical and Euclidean wavelets states that the inverse stereographic 
projection of an admissible wavelet on the plane yields an admissible wavelet on the sphere. Mother spherical 
wavelets may therefore be constructed from the projection of mother Euclidean wavelets defined on the plane: 



where ?/^^2 G L^(R^) is an admissible wavelet in the plane. Directional spherical wavelets may be naturally 
constructed in this setting - they are simply the projection of directional Euclidean planar wavelets on to the 
sphere. In Fig. [3] we illustrate the three spherical wavelets that have been used to detect the ISW effect: 
the spherical Mexican hat wavelet (SMHW); the spherical butterfly wavelet^ (SBW); and the spherical second 
Gaussian derivative wavelet (S2GDW). Each spherical wavelet is constructed by the stereographic projection 
of the corresponding Euclidean wavelet on to the sphere. The SMHW is proportional to the Laplacian of a 
Gaussian, the SBW is proportional to the first partial derivative of a Gaussian in one direction and the S2GDW 
is proportional to the second partial derivative of a Gaussian in one direction. See Ref. 25 for formal definitions 
of these wavelets. 

3.3 Steerability 

The wavelet transform outlined above is inherently directional (since it lives in the rotation group SO (3)), allowing 
for the construction of both axisymmetric and directional wavelets (an axisymmetric wavelet is invariant under 
rotations around itself; all wavelets that are not axisymmetric, we refer to as directional). For example, the 
SMHW (Fig. [3^) is axisymmetric, whereas the SBW (Fig.[3lD) and the S2GDW (Fig.[3t) are directional. Certain 
classes of directional wavelets are said to be steerable. A steerable wavelet for any arbitrary rotation about itself 
by 7 may be represented by a weighted linear combination of non-rotated basis wavelets. In some cases the basis 
wavelets are simply rotated versions of the steerable wavelet itself. Due to the linearity of the wavelet transform 

^Sometimes referred to as the spherical first Gaussian derivative wavelet. 
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(a) SMHW (b) SBW (c) S2GDW 

Figure 3. Spherical wavelets at scale a = 6 = 0.2. 



this property may be transferred to the wavelet coefficients. For steer able wavelets it is possible to compute the 
wavelet coefficients for any continuous orientation 7 from a weighted sum of wavelet coefficients computed using 
the non-rotated basis wavelets: 



where t/jm are the M basis wavelets and kmij) are the interpolating functions. Both the SBW and the S2GDW 
are steer able, requiring M = 2 and M = 3 basis wavelets respectively. For a more in-depth discussion of steer able 
wavelets on the sphere, with specific examples and interpolating functions specified, see Refs. 22,25. 

The ability of steerable wavelets to probe all continuous orientations mean that it is possible to extract the 
single orientation corresponding to the most dominant feature at each position on the sphere. It addition to 
the orientation of this feature, it is also possible to compute other morphological measures of the feature. For 
example, the S2GDW may also be used to compute measures of the signed-intensity and elongation of local 
features. We denote the morphological measures computed at each position on the sphere as /, D and E for 
the signed-intensity, orientation and elongation respectively. By allowing one to probe the morphology of local 
features, steerability facilitates an additional approach that may be used to probe dark energy. 



If dark energy exists, a correlation between the CMB and the LSS of the Universe is predicted by the ISW effect. 
In the absence of dark energy no correlation should exist. To verify the existence of dark energy one simply needs 
to detect a correlation between the CMB and a tracer of the LSS, such as the nearby radio galaxy distribution 
(and rule out other sources for any observed correlation). Wavelets are used to search for this correlation due to 
the ability to simultaneously localise signal components in scale and position, thereby improving the sensitivity 
of the analysis. In the remainder of this section we discuss the analyses performed in Refs. 12, 14, 15 and present 
the results obtained. Firstly, we describe the methods used, detailing the wavelet estimators applied, the data 
examined and the analysis procedures followed. We then present the detections of dark energy that are made 
and the constraints that are placed on the properties of dark energy. 

4.1 Wavelet estimators 

The correlation between the CMB and LSS data is computed in wavelet space by 
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4. DETECTING AND CONSTRAINING DARK ENERGY 
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where Si = {W^ I ^ E} is the term correlated (wavelet coefficient or morphological signed-intensity, orientation 
or elongation respectively), with superscript N or T representing respectively the LSS and CMB data and the 



(a) WMAP (b) NVSS 

Figure 4. WMAP co-added three- year and NVSS maps after appUcation of the joint mask. 



bar denotes a mean computed over the sky (note that the correlated terms are computed at A^p discrete points 
only). The scale parameter is given by a = (a, 6)^, which is represented by a only when an isotropic analysis 
is performed. In the absence of an ISW effect the CMB and LSS data should be independent and none of the 
correlation estimators defined by Eqn. ([16]) should exhibit a significant deviation from zero. It is possible to 
define estimators of the correlation in real and harmonic space also. However, it has been shown that the wavelet 
estimator is optimal for a large range of scales. The predicted effectiveness of various wavelets for detecting 
the ISW effect has also been compared, concluding that the wavelets that have been adopted are expected to 
perform well.^^ 

4.2 Data and simulations 

The three- year release of the WMAP data and the NVSS radio galaxy survey are examined for correlations 
induced by the ISW effect. WMAP contains many receivers with different channels and frequency bands in order 
to eliminate foreground contamination and instrumental systematics from observations. In order to enhance the 
signal-to-noise ratio, a co-added WMAP map is constructed from a noise weighted sum of the maps observed 
by different channels and receivers. Residual Galactic emission and bright point sources are removed by the 
application of a mask. It is this co-added WMAP map that is examined for correlation with the NVSS data. 
The large sky coverage and source distribution of the NVSS data make it an ideal probe of the local matter 
distribution to use when searching for the ISW effect. Sources in the catalogue are thought to be distributed in 
the redshift range < z < 2, with a peak distribution at z ~ 0.8.^ This corresponds closely to redshift regions 
where the ISW signal is expected to be produced. Not all of the sky is sufficiently observed in the NVSS 
catalogue. A joint mask is constructed to exclude from subsequent analyses those regions of the sky not observed 
in the catalogue, in addition to the regions excluded by the WMAP mask. The WMAP and NVSS data, with 
the joint mask applied, are illustrated in Fig. [H 

Monte Carlo simulations are performed in order to constrain the statistical significance of any detection of 
correlation between the WMAP and NVSS data. 1000 simulations of the WMAP co-added map are constructed, 
modelling carefully the beam and anisotropic noise properties of each of the WMAP channels and mimicking the 
co-added map construction procedure. No correlation exists between these WMAP simulations and the NVSS 
data, hence the analyses subsequently performed on the data may be repeated on the simulations and compared 
in order to constrain the significance of any correlation detected in the data. 

4.3 Procedures 

The analysis procedures consist of computing the wavelet correlation estimators described by Eqn. (p!6|) for the 
various cases. In the analysis performed by Ref. 12 the correlation is computed from the wavelet coefficients of the 
co-added WMAP and NVSS data using the axisymmetric SMHW for a range of isotropic scales. The correlation 
is computed from the wavelets coefficients of both the SMHW and the SBW by Ref. 14. In this directional setting 
anisotropic dilations are adopted and a range of 7 orientations are considered. The steerable S2GDW is applied 
in Ref. 15 to compute morphological measures of local features, which are then correlated. Any deviation from 
zero in any of these estimators is an indication of a correlation between the WMAP and NVSS data and hence 
a possible detection of the ISW effect. An identical analysis is performed using the simulated co- added maps in 
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Figure 5. Wavelet correlation surfaces. Contours are shown for levels of two and three A^^- 

place of the WMAP data in order to construct significance measures for any detections made. Finally, for the 
analyses that do not examine the morphology of local features, it is possible to use detections of the ISW effect 
to constrain properties of dark energy. 

4.4 Results 

Detections of the ISW effect and, consequently, evidence for the existence of dark energy are presented in this 
section. We first concentrate on the wavelet coefficient correlations performed by Refs. 12,14, before turning our 
attention to the correlation of measures of the morphology of local features, facilitated by the steerable S2GDW, 
performed by Ref. 15. 

4.4.1 Wavelet coefficient correlation 

Highly statistically significant detections of the ISW effect have been made using the wavelet correlation approach 
first pioneered by Ref. 12, and since extended to a directional analysis by Ref. 14. The distribution of the wavelet 
correlation statistics were found to be Gaussian distributed, hence the approximate significance of detections of 
correlation may be inferred directly from the number of standard deviations (denoted Ncr) that the correlation 
observed in the data deviations from the mean value computed from Monte Carlo simulations. Ncr surfaces 
computed for the SMHW and the SBW are shown in Fig. [5] in the anisotropic dilation space. For both wavelets, 
the maximum correlation detected is made at a level of N(j — 3.9 on wavelet scales about a = (100, 300)^ 
arcminutes.^^ It should be noted that the level quoted for these detections is based on an a posteriori selection 
of the scale corresponding to the most significant detection. Quantifying the significance in this manner is 
satisfactory, provided that the a posteriori nature of the analysis is acknowledged. The sign (positive) and the 
scale of the detected correlation are consistent with an ISW signal. Moreover, foreground and instrumental 
systematics have been analysed in detail and were found not to be the responsible for the correlation detected 
with either wavelet. As a final consistency check, the regions on the sky contributing most significantly to 
the correlation observed were localised and examined in closer detail. These regions were determined not to be 
atypical and do not, in general, correspond to regions with particularly bright point sources that potentially 
could have contaminated the data.^^ These tests strongly suggest that the correlations detected by Refs. 12,14 
are due to the ISW effect, thereby independently verifying the existence of dark energy. 

With the ISW effect detected successfully, it is possible to use the detection to constrain cosmological param- 
eter that describe dark energy, such as the proportional energy density Qa and the equation of state parameter 
w. Current constraints on suggest that dark energy constitutes approximately 74% of the energy density of 
the Universe, i.e. Qa — 0.74, while constrains on w are consistent with a value of —1 (corresponding to the simple 
cosmological constant case).^ By comparing the wavelet correlation measured from the data with theoretical 
prediction^ (which vary as a function of Qa and w)^ it is possible to compute the likelihood distribution for the 



"The theoretical wavelet correlation for isotropic dilations is derived by Ref. 12, while the extension to anisotropic 
dilations is derived by Ref. 14. 
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Figure 6. Likelihoods constructed using the SMHW for parameters {Qa,w). The full likelihood surface is shown in 
panel (a), with 68%, 95% and 99% confidence contours also shown. Marginalised distributions for each parameter 
are shown in the remaining panels, with 68% (yellow/light-grey), 95% (light-blue/grey) and 99% (dark-blue/dark-grey) 
confidence regions also shown. The parameter estimates made from the mean of the marginalised distribution are shown 
by the triangle and dashed line. 

parameters. In the case of uniform priors on the parameters, the likelihood is identical to the posterior distri- 
bution, and parameter estimates can be inferred from the likelihood directly. In Fig. [6] the likelihood recovered 
using the SMHW is illustrated (the likelihood recovered using the SBW differs slightly, although inferences drawn 
from it are consistent with the SMHW case). Marginalised distributions for Qa and w are recovered from the 
joint likelihood, with parameter estimates of Qa = 0.63^o 17 

and w = -0.77t°oil computed from the mean of the 
marginalised distributions. These parameter estimates are consistent with estimates made by numerous other 
analysis techniques and data sets. Although wavelets perform very well when attempting to detect the ISW effect 
since one may probe different scales and positions, once all information is incorporated to compute the likelihood 
surface the performance of a wavelet analysis is comparable to other linear techniques, as expected. The power 
of the wavelet analysis is the ability to make highly significant detections in the first instance. Nevertheless, the 
constraints placed on dark energy parameters provide a coherent picture, further verifying the existence of dark 
energy. 

4.4.2 Local morphological correlation 

A local morphological analysis has been performed by Ref. 15 recently to search for the ISW effect, resulting in 
an alternative highly significant detection. Rather than correlating the wavelet coefficients themselves, measures 
of the morphology of local features are correlated in this analysis, as described previously. Not only does this 
probe the existence of the ISW effect, and hence dark energy, but it also provides insight into the morphological 
nature of the correlation induced between the CMB and LSS. This further insight may in future help to better 
understand the nature of dark energy. The correlation of morphological measures of signed-intensity, orientation 
and elongation computed from the WMAP and NVSS data are displayed in Fig. [7j with significance levels 
obtained from Monte Carlo simulations also shown. Only isotropic dilations are considered in this analysis. A 
significant detection of correlation is made in the signed-intensity of local features, with the correlation computed 
from the data falling outside of the 99% significance level on a scale of a = 400 arcminutes. Moderate detections 
are also made in the other morphological measures (orientation and elongation), although the significance of 
these detections is lower. In this analysis a test was performed to examine the significance of the correlation 
when the statistics for all scales are considered in aggregate (as opposed to the a posteriori analysis performed 
when considering a single scale only). Using this test, correlation is detected in the signed-intensity of local 
features at 95% significance. On the scale (a = 400 arcminutes) corresponding to the most significant detection 
of correlation in the signed-intensity of local features, it may even be possible to observe the correlation by eye 
(see Fig. [8j). Foregrounds and instrumental systematics were again examined and were determined not to be 
the source of the correlation signal observed. The correlation detected indeed appears to be due to the ISW 
effect, thereby providing an alternative detection of dark energy. In this setting, the use of the detection has not 
yet been used to constrain dark energy properties. The analysis is complicated by allowing the orientation of 
wavelet coefficients to vary as a function of the data, hence the derivation of the theoretical correlation for each 
local morphological measure is not easily tractable. The derivation of theoretical correlations and dark energy 
constrains in this setting is the focus of current research. 
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(a) Signed- intensity (b) Orientation (c) Elongation 

Figure 7. Correlation statistics computed for each morphological measure. Significance levels obtained from 1000 Monte 
Carlo simulations are shown by the shaded regions for 68% (yellow/light-grey), 95% (magenta/grey) and 99% (red/dark- 
grey) levels. 




(a) WMAP (b) NVSS 

Figure 8. Morphological signed-intensity maps corresponding to the scale on which the maximum detection of correlation 
is made. Arguably, it may be possible to observe some correlation between maps (a) and (b) by eye. 

5. CONCLUSIONS 

The ISW effect induces a correlation between the anisotropics of the CMB and the LSS of the Universe that 
exists only in the presence of dark energy. This provides an independent physical phenomenon that may be used 
to verify the existence of dark energy. Moreover, the particular correlation signal induced by the ISW effect 
may, in certain cases, be used to constrain properties of dark energy. The use of wavelets on the sphere to 
detect a correlation between the WMAP CMB data and the NVSS radio galaxy survey has been presented. Two 
approaches have been reviewed. The first approach searches for a correlation between the wavelet coefficients of 
two data sets,^^'^^ whereas the second searches for a correlation between measures of the morphology of local 
features in the two data sets.^^ Positive detections of the ISW effect are made using both approaches, thereby 
independently verifying the existence of dark energy. The effectiveness of a wavelet analysis on the sphere is 
demonstrated by the high statistical significant of the detections. Various tests have shown that in all cases the 
correlation detected is indeed due to the ISW effect and not contamination in the data or instrumental effects. 
Furthermore, for the first of the two techniques (the wavelet coefficient correlation), it is possible to use the 
correlation signal detected to constrain cosmological parameters that describe dark energy (using the correlation 
detected in the morphological analysis to constrain dark energy parameters is the focus of current research). 
Parameter estimates have been obtained that are consistent with estimates made using numerous other analysis 
techniques and data sets. The independent analyses performed therefore provide a coherent picture, further 
verifying the existence of dark energy. 

In more sophisticated dark energy models, perturbations in the dark energy fluid should be taken into account. 
Wavelet based ISW detections could be used to probe this property and to constrain the corresponding sound 
speed of dark energy (see e.g. Ref. 27). This is an interesting area of future research, where again the use of 
wavelets on the sphere to detect the ISW effect may help us to better understand the nature of dark energy, the 
dominant component of our Universe of which we know very little currently. 
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